Density ripples in expanding low-dimensional gases as a probe of correlations 
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We investigate theoretically the evolution of the two-point density correlation function of a low- 
dimensional ultracold Bose gas after release from a tight transverse confinement. In the course 
of expansion thermal and quantum fluctuations present in the trapped systems transform into 
density fluctuations. For the case of free ballistic expansion relevant to current experiments, we 
present simple analytical relations between the spectrum of "density ripples" and the correlation 
functions of the original confined systems. We analyze several physical regimes, including weakly 
and strongly interacting one-dimensional (ID) Bose gases and two-dimensional (2D) Bose gases 
below the Berezinskii-Kosterlitz-Thouless (BKT) transition. For weakly interacting ID Bose gases, 
we obtain an explicit analytical expression for the spectrum of density ripples which can be used 
for thermometry. For 2D Bose gases below the BKT transition, we show that for sufficiently long 
expansion times the spectrum of the density ripples has a self-similar shape controlled only by the 
exponent of the first-order correlation function. This exponent can be extracted by analyzing the 
evolution of the spectrum of density ripples as a function of the expansion time. 

PACS numbers: 03.75.Hh, 67.85.-d 



I. INTRODUCTION 

A. Quantum noise studies of ultracold atoms 

Quantum correlations can be used to identify and 
study interesting quantum phases and regimes in ultra- 
cold atomic systems. Recent experimental advances in- 
clude detection of the Mott insulator phase of bosonic [1] 
and fermionic [H atoms in optical lattices, production 
of correlated atom pairs in spontaneous four-wave mix- 
ing of two colliding Bose-Einstein condensates 3], stud- 
ies of dephasing [4] and interference distribution func- 
tions |5(] in coherently split one-dimensional (ID) atomic 
quasicondensates (QC), observation of the Berezinskii- 
Kosterlitz-Thouless (BKT) transition [1, Q in two- 
dimensional (2D) quasicondensates Q, and Hanbury- 
Brown-Twiss correlation measurements for nondegen- 
erate (ND) metastable "^He JS] and ^He atoms [10| . 
bosonic and fermionic |12l | atoms in optical lat- 
tices, and in atom lasers 11 3l . In one-dimensional atomic 
gases [l3, [l^ [H, [13, [H, Il9l.[20j , in situ measurements of 
correlations have been attained by means of photoasso- 
ciation spectroscopy [2]| or by measuring the three-body 
inelastic decay Q , using the proportionality of the cor- 
responding rates to the zero-distance two-particle and 
three-particle correlation functions, respectively [23j . 

Recently it was demonstrated that one can detect sin- 



gle neutral atoms in a tight trap or guide [2J, |25|, [26!, |27|, 
[28l . [29j . However, direct (not inferred from any kind of 
atomic loss rate [2l|, [2^ ) observation of interatomic cor- 
relations at short distances in trapped ultracold atomic 
gases is hindered in many cases by either the finite spa- 
tial resolution of the optical detection technique or the 
very low detection efficiency of the scanning electron mi- 
croscope [27|. Therefore one needs to release ultracold 
atoms from the trap, diluting the atomic cloud in the 
course of expansion. 

In this paper we address the question of how the cor- 
relations in the low-dimensional system evolve during 
the time-of-flight expansion, and discuss how the den- 
sity variations in the time-of-flight images relate to the 
properties of the original trapped quantum gas. These 
"density ripples" in the expanding gas reflect the origi- 
nal thermal or quantum phase fluctuations existing in the 
cloud under confinement. Such phase fluctuations are al- 
ready present in three-dimensional (3D) Bose-condensed 
clouds under an external confinement with large aspect 
ratio [1^. Their effect on density ripples of expanding 
clouds has been observed [3l|, |32| . [33| , but quantitative 
analysis of such experiments was complicated since one 
had to take into account interactions in the course of 
expansion. However, for sufficiently strong transverse 
confinement reached in current experiments with low- 
dimensional gases (chemical potential of the order of 
the transverse confinement frequency), the gas expands 
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rapidly in the transverse direction so interactions dur- 
ing the expansion stage can be safely neglected. Then 
one can develop a simple analytical theory, which di- 
rectly relates the spectrum of the density ripples after 
the expansion to the correlation functions of the origi- 
nal fluctuating condensates. Similar question has been 
considered for 3D clouds expanding in the gravitational 
field but only for noninteracting atoms [s^l . We also note 
the density ripples we discuss are different from the den- 
sity modulations which appear due to interactions during 
expansion and have been studied in Refs. [s^ and [36j . 



B. Density ripples in expanding condensates: 
Preview 

We consider one- or two-dimensional atomic gases re- 
leased from a tight trap formed by a scalar potential 
as realized on atom chips or in optical lattice experi- 
ments. We consider the situation when free expansion 
takes place in all three dimensions. This should be con- 
trasted to the expansion of such a gas inside a waveguide 
[H M, m, m, So, El, m m Si, 0] , with the transverse 
confinement being permanently maintained. In the latter 
case, the nonlinear atomic coupling constant 
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where lo± is the transverse trapping frequency and is 
the atomic s-wave scattering length, remains the same. 
While a bosonic gas rarities during such expansion, col- 
lisions remain important. For example, in the ID case 
dynamics asymptotically reaches the limiting Tonks- 
Girardeau (TG) [i^ regime of impenetrable bosons. In 
our case, if the fundamental frequency of the potential 
of the transverse confinement is much larger than the 
initial chemical potential of the atoms, the expansion in 
the transverse directions is determined mainly by the ki- 
netic energy stored in the initial localized state of the 
transverse motion. Interatomic collisions play almost no 
role in the expansion. Moreover tight transverse con- 
finement decouples the motion of trapped atoms in the 
longitudinal and transverse directions. Thus when ana- 
lyzing density ripples we can reduce the problem to the 
same number of dimensions as the initial trap (see dis- 
cussion below in Sec. |TT]). For a ID trap we consider 
a one-dimensional spectrum of density ripples, and for 
atoms which were originally confined in a pancake trap 
we analyze two-dimensional density ripples. 

Before we consider a general formalism, it is useful 
to present the analysis for the simplest situation. Let 
us assume that the initial state can be described using 
the mean- field Bogoliubov approach f47, 48, 49]. Let 
be the creation operator of atoms at momentum k 

k 

right before the expansion. After free expansion during 
time t, in the Heisenberg representation we have ^^(0 — 

"^t e* 2m ^ where m is the atomic mass. Then the density 



operator at time t is given by 



(2) 



and for the density correlation function we obtain 
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The expectation value ) should be taken 

f ^ fei fe2 fea kil 

in the original condensate before the expansion. Within 
the mean-field Bogoliubov theory only a state with k = Q 
is macroscopically occupied. Thus in Eq. ([3]) we take 
two operators to be ^/N = y/nujL, where tlid is the 
atomic density before the expansion. Thus Eq. ([3]) can 
be written as 
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The Bogoliubov theory predicts expectation values of 
(*%*t), and 1 + 2(*t*,) as 
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where eq — h q /{2m), fj, = gionio is the chemical po- 
tential, Eq — -y/ eq{2^ + eq) is the Bogoliubov excitation 
spectrum, and Hb is the Bose occupation number. 

From these equations we can easily find the mean-field 
spectrum of density ripples {\p^^^ {q)\'^) [see Eq. (PT|) and 
the discussion nearby for the precise mathematical defi- 
nition of the spectrum] 
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The general character of the spectrum is clear from 
Eq. ([7|). As a function of momentum it is not mono- 
tonic. We find minima near hq'^t/m — 27rn, and maxima 
close to hq-t/m — 7r(2n — 1), where n is a positive in- 
teger number. Note that while the positions of maxima 
and minima are essentially universal, the amplitude of 
individual maxima depends on both temperature and in- 
teraction strength. 

The mean-field analysis leading to Eq. ([7]) is concep- 
tually simple but has limited applicability. It is appli- 
cable for weakly interacting 3D Bose condensates if the 
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interactions during expansion are switched off using Fes- 
hbach resonances i5Q] . In lower dimensions, thermal and 
quantum phase fluctuations are expected to suppress true 
long range order in 2D Bose condensates at finite tem- 
perature |i5l| and in ID Bose condensates even at zero 
temperature [s^ ]. In this paper we show how the analy- 
sis of the density ripples can be extended to more com- 
plicated but experimentally relevant situations when the 
mean-field approach breaks down. We will find a simi- 
lar structure to Eq. ([7]): positions of maxima and min- 
ima of the spectrum are given by the same approximate 
universal conditions on the momenta. However explicit 
expressions for the strength of individual maxima will be 
very different. They will contain rich information about 
fiuctuations of low-dimensional condensates. 



C. Relation to other work 

Conceptually, the question we consider in this paper 
is somewhat similar to the interpretation of cosmolog- 
ical observations. In the latter case, quantum fluctu- 
ations present in the early universe after its expansion 
result in observable anisotropics of the cosmic microwave 
background radiation ^53; , .54] , and in the density ripples 
of matter which eventually evolve into galaxies |55j . In 
our case, density ripples of the expanding clouds contain 
important information about correlations present in the 
trapped state. Analogies between properties of conden- 
sates and cosmology have attracted significant attention 
recently [H [57, 58, 5^ M, IMI, [13, [H . 

In addition to the mentioned approaches, several other 
techniques have been used to experimentally study cor- 
relations of low-dimensional gases. Some of them rely on 
creation of two copies of the same cloud [6^ [gsI . \qK 67 , 
[ssj , while others require analysis of noise correlations 69 1 
or in situ density-fluctuation statistics [tO]. Interference 
experiments between two low-dimensional clouds [1, Q 
can also be used to characterize two- poin t and multi- 
point correlation functions [tH. [t^. [tsI. |74|. Analysis of 
density ripples is a much simpler experiment and, as we 
discuss in this paper, can be used for thermometry. This 
is particularly important for weakly interacting ID Bose 
quasicondensates [tI, I?^, for which the standard ap- 
proach to measuring temperature by fltting density pro- 
files cannot be extended to temperatures of the order of 
the chemical potential. In this regime, the chemical po- 
tential is very weakly dependent on the temperature [73| , 
thus finite temperature leads only to small corrections 
to the "inverted parabola" density profile [78J. An im- 
proved thermometry method based on comparison of in 
situ measured density profiles with solutions of Yang and 
Yang equations [t^ in the local-density approximation 
has been developed in Ref. fl^ . 

There has been significant theoretical interest in cor- 
relation functions of the ID Bose gas. At distances much 
larger than the healing length, correlation functions are 
described by Luttinger liquid theory [13, [8l|, [111 . In the 



weakly interacting quasicondensate regime, correlation 
functions can be described by extension of Bo goliu bov 
theory to low-dimensional gases [zl, [tI [H, [Mill- In 
the strongly interacting regime, one can use "fermion- 
ization" [46] of a ID Bose gas to evaluate correlation 
functions at all distances as certain determinants [s^ . 
The Lieb-Liniger model [13] which describes the ID Bose 
gas is exactly solvable, and one can also analytically ob- 
tain zero-distance two-point [H, [s^ and three-point [13] 
density correlations for any interaction strength, and ex- 
tract certain dynamical correlation functions [9l|, [H, [H, 
[9^ . [9^ from the exact solution. Various numerical tech- 
niques have been used as well [zl, [H, [13, [11] , and recent 
results including the decoherent quantum regime [H, [9§] 
are summarized in Refs. 
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D. Structure of the paper 

This paper is organized as follows. In Sec. [Tlj we de- 
rive simple analytical relations between the density rip- 
ples after the expansion and the correlation functions of 
the original system before the expansion. In Sec. IIII Al 
we analyze the case of weakly interacting ID Bose gases 
and obtain explicit expression for the spectrum of density 
ripples. In Sec. IIIIBI we consider the case of a strongly 
interacting ID Bose gas. In Sec. IIII CI we review general 
features of the density-density correlation function in ex- 
panding ID Bose clouds. In Sec. |TV] we discuss 2D Bose 
systems below the BKT transition [6[, [3] . We summarize 
our results and make concluding remarks in Sec. [Vj 



II. FREE EXPANSION 

In this section we focus on the atoms expanding from a 
onc-dimensional trap. The atom field oper ator evolution 
during the free expansion is given by [ll5l ] 



*(r,0 = / dVG3(r-r',t)*(r',0), 



(8) 

where the Green's function of free motion is 
G3(r-r',t) = Gi{x-x',t)G,{y-y',t)G,{z-z',t), (9) 



2iTi ht 



exp I 



2U 



(10) 



with m being the atomic mass. Tight transverse con- 
finement decouples the motion of trapped atoms in the 
(y, z)-plane and along the waveguide axis x, so that the 
transverse motion is confined to its ground state z), 
and ^'(r,0) = f±{y,z)ip{x,0). This, alongside with Eq. 



4 



([9]) , allows for a separation of motion in the longitudinal 
and transverse directions, effectively reducing the prob- 
lem to ID. 

We introduce the two-particle density matrix for the 
longitudinal motion as 

p{xi,X2;x[,x'2\t) = I^%i)\x[,t)%i^\x'2,t)-li}{x2,t)'lj}{xi,ty^ 

(11) 

Then we define the two-point density correlation function 

p{Xi,X2\Xi,X2;t) 



where 
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where n{x,t) = yijj^ {x,t)ili{x,t)J . The free evolution of 

the two-particle density matrix is given by the convolu- 
tion of the two-particle density matrix at i = with four 
respective Green's functions, one for each spatial argu- 
ment, two of the Green's functions being complex conju- 
gate. We are interested in the case xi = x'l and X2 — x'2 
in the final state. Then we obtain 

p{xi,X2;xi,X2:,t) ^ 

dxz / dx'^ I dx4 / dx'4^Gi{xi — X3,t)Gi{x2 — X4,t) 



xG'Kxi - x'3,t)Gl{x2 ~ x\,t)p{xzTX4\x'n,,x\;Q). (13) 

We assume that the product of the typical velocity of 
the atoms in the a;-direction and the expansion time is 
much smaller than the size of the trapped atomic cloud. 
Then we are allowed to consider a uniform sample with 
length L 00, with the ID number density nuj = N/L 
being kept constant in the thermodynamic limit {N being 
the total number of atoms). Note that this limit is oppo- 
site to the conventional limit of infinitely large expansion 
times, in which density in real space reflects the initial 
momentum distribution (in that regime, it was recently 
proposed [116] that noise correlations in density pro- 
files can be used to probe properties of low-dimensional 
gases). 

In our limit n{x, t) — nio is constant in time, and 
the two-particle density matrix is translationally invari- 
ant (it does not change if all four of its spatial arguments 
are shifted by the same amount) at any time. The den- 
sity correlation function then depends on the coordinate 
difference only, so we use the notation 



52(2:1 - X2;t) = g2{xi,x2;t). 



(14) 



Using the translational invariance of the two-particle den- 
sity matrix and the identity 
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we arrive at 
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-{xi-X2-x')^]}g{x;x';0), (16) 
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Obviously, p{xi,X2;xi,X2]t) = q{xi — X2]Xi — X2]t) = 
"ii532(a;i — X2]t). The physical meaning of Eq. (fT6|) is 
that the motion of the center of mass of an atomic pair 
plays no role in the dynamics of establishing g2{xi—X2] t), 
which is fully determined by the relative motion. The 
relative-motion deg ree o f freedom is characterized by the 
reduced mass m/2 (ll7j . 

Let us now consider some properties of the two-particle 
density matrix g{xi; X2] t) for bosons. Changing the sign 
of xi or X2 is equivalent to a permutation of two bosons 
and, hence, does not change the two-particle density ma- 
trix, i.e.. 



g{xi\x2\t) = g{\xi\\ \x2\]t). 



(18) 



For the regimes we consider the density matrix of neutral 
bosons can be assumed to be real. This, together with 
the Hermicity property, results in 



g{xi]X2;t) = g{x2]xi;t). 



(19) 



Using Eqs. IjlSp and (|19p and Fourier transforming the 
Green's functions, we can reduce Eq. to 



g{x; x; t) 



00 poo 



dq / dX cos qx cos qX 
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xgi\X-^\-\X + ^\;0). (20) 
m m 



Alternatively, this equation can be written as 
{\p{q)\')^C^dXcosqX 
X /v;t(M, o)^t(x, 0)V3(X + 0)^(0, 0)\ . (21) 



Here (|p(g)P) is the spectrum of density ripples at time 
t, which in experiment can be obtained by Fo urier trans- 
forming absorption images after expansion [ll8j |. It is 
related to two-point density correlation function as |ll9l | 



{\p{q)?) ^ rilj, dxe,^^{iqx){g2{x-t)-l). (22) 

J —00 

Equations (PO)) and (PT|) provide a simple analytical re- 
lation between the properties of the density ripples after 
the expansion and the correlation functions before the 
expansion. 

It is straightforward to generalize the above analysis 
to the 2D case. In particular, the analog of Eq. (PTjl for 
the time evolution of the two-point density correlation 
function has the same form, with X substituted by r and 
q treated as a 2D vector. Namely, we obtain 



(|p(q)P) = Sr2 rf^rcosqr 
X (V'^^, 0)^-^1-, 0)^(r + ^,0)^(0, 0)) 



(23) 
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III. ID BOSE GASES 
A. Weakly interacting ID Bose gases 

In this subsection we will consider the spectrum of den- 
sity ripples of weakly interacting ID quasicondensates. 
Before we proceed to full analytical theory, let us return 
to the simple mean-field Bogoliubov approach, which we 
discussed in the introduction. Readers may be skeptical 
about the applicability of the mean-field approach to ID. 
Indeed, there is no long-range order for ID gases even at 
zero temperature [52|, and mean-field approach is gen- 
erally not applicable. However, in the regime of weak 
interactions and under certain conditions on the expan- 
sion time t, the spectrum of density ripples is captured 
correctly by the mean-field approach, as we will verify 
later in a rigorous calculation. 

We consider Eq. (O in the hmit <C fi,Eq ^ ksT. 
In this case one can neglect the first term in the second 
parentheses, use an approximation Eq cx g, and expand 
the Bose occupation number, leading to 



no 



2mkBT [1 - cos 



a — 



(24) 



Let us now present a full calculation, which does not 
make a mean-field approximation. For weak interactions 
Bogoliubov theory has been extended to low-dimensional 
quasicondensates [t^, and can be used to calculate cor- 
relation functions at all distances. For quasicondensates, 
the fiuctuations of the phase are described by the Gaus- 
sian action. For Gaussian actions, higher order correla- 
tion functions are simply r elate d to t wo-p oint correlation 
functions (see, e.g., Refs. [l20j | and [l2l( ). and the four- 
point correlation function in Eq. (I2ip factorizes into prod- 
ucts of two-point correlation functions of bosonic fields 
as fzl 



V't(a;'i,0)V^(4,0)V(a;i,0)V(x2,0) 



^t(a;^,0)^(x„0) 



V^t (x'l , OMx^ , 0) ) (v-t (a;i , 0)^(0.2 , 0) 



(25) 



This equation gives correct result for all values of 
xi,X2tx'i, and x'2 in the leading order over \/K expan- 
sion (see the definition of _ftr ^ 1 below), and can be used 
to evaluate the spectrum of the density ripples in weakly 
interacting condensates for all times. The two-point cor- 
relation function ^'(/'^(x']^, 0)'!/'(a;i, 0)^ = niugi{x'i~xi) is 

translationally invariant, and is simply related to predic- 
tions of the Bogoliubov theory. For ID quasicondensates, 
one has [ttI 



^t(x, 0)^(0,0) 



niD 



exp 



2K^ \Ch 
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where K = Trhnnj / (mc) ^ 1 is the Luttinger liquid pa- 
rameter, c is the speed of sound, = h/^^rnjl is the heal- 
ing length, and /i is the chemical potential. When only 
lowest transverse mode is occupied (/i <C hoj±), speed 
of sound is given by c = yj2Tujj±niDas/ra. The dimen- 
sionless function /(s) depends on the temperature, and 
equals 

/(s) = 2 / dk{l-cosks){[ul + vl]nk + vl], (27) 

JO 



where 



Uk 




(28) 

(29) 
(30) 



For finite temperatures, the function /(s) has the follow- 
ing asymptotic behavior 

/(5) ^ 7r|s|-^ + C for 7r|s|^ » 1, (31) 

where C = C{kBT/^) is of order 0(1) for knT^^ fj,. 

Quasicondensate theory is valid [73, llOOL Il0l| for tem- 
peratures 



keT/n < K/tt, 



(32) 



significantly beyond the regime of validity of Luttinger 
liquid theory, which is restricted to kgT/fi <^ 1. The 
longitudinal density profile of a quasicondensate in exter- 
nal harmonic confinement follows the inverted parabola 
shape under condition ([5^ . see, e.g., Ref. [zl]. Due to the 
low fraction of the thermally populated excited states, it 
is problematic to extract the temperature of the gas from 
fitting bimodal distributions to the observed density pro- 
files. Below we show that the spectrum of density ripples 
can be used as a convenient tool to characterize the tem- 
perature, and is sensitive to temperatures of the order of 
the chemical potential. 

To be specific, let us consider the case of ^'''Rb 
atoms (scattering length — 5.2 nm) with density 
riiD — 40 /im~^ and transverse confinement frequency 
Lu± = 2tt X 2 kHz, resulting in Luttinger liquid parame- 
ter X sa 47 and healing length S^h ~ 0.37 /im. We can use 
Eqs. to numerically evaluate in-trap correlation 

functions. By performing then a numerical integration of 
Eq. (|2ip for various temperatures and expansion times, 
we can evaluate the spectrum of density ripples under 
condition ([5^ . and the results are shown in Figs. [T]and 
[21 In the inset to Fig.[2]we also show g2{x; t) evaluated us- 
ing the inverse of Eq. . In the quasicondensate regime 
the behavior of 52(2;; t) follows the qualitative discussion 
of Sec. [mCl 
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FIG. 1: Normalized spectrum of density ripples 
{\pil)\^} /i^iD^h) for weakly interacting ID quasicon- 
densate of *'^Rb atoms with density nm = 40/im~^, 
transverse confinement frequency uj± — 2tt x 2 kHz, Lut- 
tinger liquid parameter K « 47, and healing length 
S,h ~ 0.37 /im. Expansion time is fixed at t = 27 ms (with 
« 11-8), and temperatures equal (top to bottom) 
T = 40nK{kBT/fi = 1), T = 27nK(fcflr//i = 0.67), 
T = 12nK (fcsT/^i = 0.3), and T = 0. Values on the axes 
of this and subsequent plots are dimensionless. Dots are 
obtained by numerical integration of Eq. (|2ip in the weakly 
interacting limit making use of Eqs. (I25|l - (|30|l . Solid lines 
correspond to analytical results [Eq. (|37|l ]. which are derived 
under condition (|33|) . 
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FIG. 2; (Color online) Normalized spectrum of density rip- 
ples (|p(g)P) /in\]j(^h) with the same parameters as in Fig.[T] 
but for a fixed temperature T — 27nK(fc_BT//i — 0.67), and 
various times of flight: t = 49ms (red, solid), t — 27ms (green, 
dashed), and t = 9.5ms (blue, dotted). Dots are obtained by 
numerical integration of Eq. (|2H) in the weakly interacting 
limit making use of Eqs. (|25p - pop . Lines correspond to an- 
alytical results [Eq. (137^ ]. which are derived under condition 
(|33p . Inset shows g2{x;t), obtained from (|p(?)P) using the 
inverse of Eq. ([22|) . 



There are several qualitative features that should be 
noted. The spectrum of density ripples is not a mono- 
tonic function, and can also have several maxima. The 
positions of the maxima only weakly depend on the tem- 
perature, and are mostly determined by the expansion 
time. The amplitude of the ripples, on the other hand, 
significantly depends both on the expansion time and 
temperature. 

Let us now derive a simple analytical expression for the 
spectrum of density ripples, which is valid in the regime 
[justified below after Eq. P^ ] 



> 1. 



(33) 



Under this condition one can use Eq. (j3H) and approxi- 
mate the two-point correlation function by 



gi{x) sa exp (-[xI/At) for > 



kBT' 



where At is defined by 



A7 



(34) 



(35) 



and does not depend on interaction strength, as long as 
Eq. (1321) is satisfied. 

Using Eqs. (12S1) and (IMD, the second hue of Eq. ^ 
can be written as 



gi(M/m)2gi(X)2 



gi{X — hqt/m)gi{X + Tiqt/m) 



-X hqt 

exp — — lor X < , 

At m 

and exp — otherwise. 



-hqt 



is responsible for 32(2^ 



The constant term exp ■ , 
00, t) = 1. Since according to Eq. we need to take 
a Fourier transform of the above expression, subtracting 
the constant on the whole interval (0, 00) does not affect 
(|p(<7)P) for q ^ 0, and we obtain 



'ID 



, ~x —hqt 

ax cos qx exp exp — - — 

At toAt 



(36) 

This integral can be evaluated in a closed form, and leads 
to an analytical answer 



{\piq)\') 



A 



-2hqt / t 2, t 2j, 

q-e^^ (Atqcos^^ -f 2sin^ 



nfoOi q(h (1 + A|,g2) 

(37) 

Note that the last equation reduces to Eq. (IM]) when 
At9 ^ 1 and hqt/mXT <C 1. Figures H] and 121 show an 
excellent agreement between the analytical result and nu- 
merical integration described earlier after Eq. (l32l) . The 
analytical result shows the same non-monotonic behavior 
as the numerical calculations. The parameter At defines 
a time scale 



tc 



6.5- 



(38) 
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after which only a single maximum persists. When sev- 
eral maxima and minima are present, their positions can 
be estimated by 



the condition 



7r(2n-l/2=Fl/2) 



(39) 



where the upper (lower) sign corresponds to the nth max- 
imum (minimum). These conditions can be understood 
as a "standing wave" conditions in Eq. and become 
more precise at lower temperatures. 

The appearance of minima and maxima in the spec- 
trum of density ripples can be understood in terms of 
matter-wave near-field diffraction. The analogous effect 
for light wav es (i n the spatial domain) is known as the 
Talbot effect [l22{ . Its matter- wave counterpart has been 
also observed in diffraction of atoms on a grating jl23l |. 
In our case, we observe near-field diffraction in the time 
domain. For each expansion time, a certain momentum 
contribution will be "imaged" onto itself, leading to a 
minimum in the spectrum of density ripples for a given 
momentum q. As compared to diffraction on a regular 
grating with a fixed period, the typical fluctuation length 
in the trapped cloud is not constant, but distributed 
around the thermal length X^- Therefore, minima in 
the spectrum appear for any sufficiently small expansion 
time, according to condition ((39)) . 

Condition (|33p can now be justified in the regime where 
(|/5(g)P) is near its largest values. In such case most of 
the contributions to Eq. (I2ip come from distances of the 
order yffit/m, and Eq. dSS]) follows from Eq. p2)) . 

So far we have been assuming that the quasiconden- 
sate is deep in the ID regime, /i, /c^T ^ Tilo^. While 
Eqs. ([25l) -([30 | are valid only under such assumption, 
Eqs. ([Ml and ([55]) also work in the weakly interacting 
quasi-lD regime, 



(40) 



Indeed, they rely only on the ID nature of long-range 
correlations, weakness of interactions, and the property 
cK = TTTiijp/m, which is a consequence of the Galilean 



mvariance |8l|. In Eqs. ([32)) and ([33)) . the Luttinger liq- 
uid parameter K can be obtained as K — hTrniD/{mc), 
where the square of the sound velocity c can be deter- 
mined from compressibility as = ni£)(9/i/9ni£))/m. 
For che mical potential /i, one can use an approximate 
relation { l24l | fj, = huj±{y/l + AasUiD — 1). 

Let us now briefly review the conditions under which 
one can neglect interactions in expanding ID clouds and 
the effects of finite condensate length L. Transverse ex- 
pansion takes place at the times of the order of inverse 
transverse confinement Wj]^. Up to the times of this or- 
der, one cannot neglect interactions during the expan- 
sion. Correlation functions which enter Eq. ([^T)) will be 
smeared up to the distances of the order Sx ~ c/uj± ~ 
^,h^J■/^±^ a-iid smearing will only weakly affect the final 
result for (|p((7)P) if qSx ^ 1. Thus to observe an oscil- 
lating spectrum of density ripples, one needs to satisfy 




« 1, 



(41) 



which easily holds for the parameters shown in Figs. [T] 
and [2 In addition, one can use Feshbach resonances [50| 
to completely switch off interactions during the expan- 
sion. 

Locally, corrections due to finite L can be neglected, 
if finite limits of integration in Eq. (|16p lead to smearing 
of delta functions up to the distances at which the cor- 
relation functions change considerably. This change can 
occur cither because of the variations of the density in 
external confinement at distances ~ L, or because of the 
decay of correlations for finite temperatures at distances 
of the order ~ K^h/a. Thus for finite temperatures these 
conditions read as 



mL 
Id 



min (i, if^h/a) > 1, 



(42) 



and are easily satisfied for parameters considered earlier, 
and, e.g., longitudinal frequency w^; = 27r x 5 Hz. Un- 
der condition ([^^ one can take the inhomogeneity of the 
density profile into account within the local-density ap- 
proximation by averaging the prediction of Eq. ([57)) . 



B. Strongly interacting ID Bose gases 

Let us now describe the evolution of the two-point den- 
sity correlation function g2{x\ t) of a strongly interacting 
ID Bose gas. A dimensionless parameter which controls 
the strength of interactions at zero temperature can be 
written as 



mgiD 2muj±as 
7 = — = ^ > 1- 



hni 



D 



(43) 



Under such conditions, the bosonic wave function takes 
on fermion properties, and the density correlation func- 
tion in the trap g2{x; 0) is the same as for non-interacting 
fermions of the same density and temperature. In par- 
ticular, it vanishes at x = 0, and one has 172(0; 0) — 0. 
However, the correlation functions that contain creation 
and annihilation operators at different points, such as 
g{xi; X2]0) in Eq. ()20[) . are not the same as for non- 
interacting fermions. This happens because bosonic op- 
erators, when written in terms of fermionic operators, 
contain a "string" which ensures proper commutation re- 
lations. 

In the Appendix we derive a representation of 
g{xi;x2',0) as a Fredholm-type determinant, which can 
be easily evaluated numerically. Combining this repre- 
sentation with Eq. (|20[) . we evaluate g2{x;t) after var- 
ious expansion times numerically. The results for zero 
temperature are shown in Fig. [31 while the results for 

finite temperature ksT — ^ ~ 1.2 *'^^2m°"* shown 
in Fig. m In spite of a considerable change in the tem- 
perature, there is no qualitative change in the behavior 
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FIG. 3: Two-point density correlation function g2{x;t) of 
a zero-temperature strongly interacting ID Bose gas (Tonks- 
Girardeau limit) for different times t after the release of the 
gas from the trap. Different curves correspond tot — (solid), 
t = 0.25m/ (hn^jj) (dashed), and t = m/{hn\]^) (dotted). 





2 3 



FIG. 4: Two-point density correlation function g2{x;t) for 
the same parameters as in Fig.[3]but for a finite temperature, 



of 52(2^; t). The qualitative behavior of g2{x; t) in Figs. [3] 
and[4]is in agreement with Eq. ((44|) below, and Ac ^ n^l) 
for the To nks- Girardeau gas. 



C. General remarks about ID case 

Before concluding this section we would like to provide 
a qualitative analysis of the evolution of the density cor- 
relation function g2{x;t) as a function of the expansion 
time t. 

The general structure of the two-particle density ma- 
trix g{xi ; 0:2; 0) of a ID Bose gas is shown schematically in 
Fig. [51 Because of the Bose symmetry, g{xi; X2', 0) is rep- 
resented in the (xi, X2)-plane by two infinite perpendic- 



xi 
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--1 
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-2 2 
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FIG. 5: Density plot of the two-particle density matrix 
£i(a;i; 12; 0) of a ID Bose gas, see Eqs. (fTTj) and (fTT)) (the 
density bar represents the q scale in units of nf^,). Initially 
(at t = 0) the bosonic system can be a Tonks-Girardeau gas 
(TG) or a weakly interacting quasicondensate (QC) or a non- 
degenerate (ND) thermal gas. The central (x\ a;2 ~ 0) part 
of the two-particle density matrix in these cases is shown in 
three respective insets. The bar shows the typical correlation 
scale Ac- 



ular "bands" of a typical transverse size Ac (correlation 
length). Asymptotically, as x\ ±00 and xi = ±2^1, 



There are several possible cases of atomic cor- 
relations near the point a;i = 12 = in a trapped ID gas. 
In general, at t = 0, we have p(0;0;0) = nf£,g2(0;0). 
In the case of the Tonks-Girardeau gas of impenetra- 
ble bosons Hi, ill 52(0; 0) EE g'^^iO) = (at zero 
temperature g2'^{x) = 1 — [sin{TTTii]jx) / {nnujx)]'^ .46:]). 
Another possibility is a weakly interacting degenerate 
gas (quasiconden sate) , where 32(0; 0) = g^'^(O) ~ 1 
[Z3, H m, [loS ^QM- Finally, the ID Bose gas can 
be non-degenerate (thermal), in which case 32(0; 0) = 
g2^ = 2. As the interparticle distance grows, the density 
correlation function quite rapidly approaches its asymp- 
totic value 32(2; ^ 00; 0) = 1 at the distances of the order 
of Ac- 

One can show that time-dependent density correlation 
function can be written as 

52(2:; t)^l + k{Xc,x, t) + [52(0; 0) - l]h{\c,x, t). (44) 

The first term (unity) stems from the band of non-zero 
values of g ahgned along the hue xi = X2 (see Fig. [5]). It 
represents the density correlation function of an ideal gas 
of distinguishable particles at equilibrium. The second 
term, K(Ac,2;,t), reflects the Bose-Einstcin statistics of 
the atoms and appears due to the second band along 
X2 = —xi. Its maximum value, k(Ac, 0, t), increases from 
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to 1 on a typical time scale ~ mX^/h. As |x| grows, this 
term asymptotically approaches on a length scale given 
by Ac- The third term describes washing-out of initial 
short-range (microscopic) correlations. The maximum 
value of h(Xc, x, t) is reached at a: = 0, it decreases from 

1 to on a time scale ^ mA^/Ti, and h{Xc,x,t) « if 
\x\ ^ Ac- In the course of free evolution, the density 
correlation properties of an expanding Bose gas become 
similar to that of an ideal Bose gas at temperature ksT ^ 



IV. 2D BOSE GASES BELOW THE 
BEREZINSKII-KOSTERLITZ-THOULESS 
TEMPERATURE 

Let us now discuss the properties of density ripples 
in expanding 2D clouds- Recently 2D condensates have 
been rea l ized experimentally in several groups 8, 15, 68, 
Il02l . Il03l Il05l | . Reduced dimensionality has dramatic ef- 
fect on thermal fluctuations- In the case of 2D Bose gases 
there is no true long-range order for any finite tempera- 
ture [in - For uniform 2D Bose clouds at sufficiently low 
temperatures, the two-point correlatio n function behaves 
at large distances as [75|, 



i^^(r,0)?/'(0,0)\«n2D 



2D 



for r > A2D. (45) 



For weakly interacting 2D Bose gas at small tempera- 
tures, one can evaluate parameters of Eq. (|45p from mi- 
croscopic theory. The dimensionless parameter charac- 
terizing weakness of interactions is written as [H, Il07l | 



87rTOaj_L 



< 1. 



(46) 



The exponent t] in Eq. (|45)) equals [85|, llQ. 



r/ = — < 1 for fcsT < fcsTd = —, 47 

Td m 

and \2D equals the de Broglie wavelength of thermal 
phonons Tic/ [kBT) at fc^T ^ and the two-dimensional 
healing length ^2_d = ?i/ sJrnpL at high temperatures 

Equation (|45p remains valid for rj smaller than 



ric = 1/4, 



(48) 



at which point the BKT [^, 0, ll2Clj | transition takes place 
due to proliferation of vortices, and correlation functions 
start to decay exponentially with distance. 

Such a transition for ul traco ld 2D Bose gases has been 
observed recently 0, [H, Il05j |. and its microscopic ori- 
gin has been elucidated. Experiments of Ref. [1] studied 
interference of two independent 2D Bose clouds, which 
requires imaging along the "in-plane" direction and in- 
evitably leads to averaging over inhomogeneous densities. 



Study of the spectrum of density ripples in expanding 
clouds with imaging in transverse direction (as done in 
Ref. [H^) avoids this problem altogether, and can provide 
access to properties of correlations at fixed density. The 
interplay between the BKT transition and the effects of 
the external confinement is a rather c omplicated question 
even for weakly interacting Bose gas llOSL 1 1041 . llOSl Il09l | , 
and we will only discuss the uniform case here. 

Even for weak interactions, one cannot use quasicon- 
densate theory to analytically describe correlations as 
functions of microscopic parameters in the vicinity of the 
BKT transition or to predict the transition te mpe rature, 
and has to resort to fully numerical methods jllOl | . Nev- 
ertheless, the factorization property [Eq. ([25]) ] remains 
valid for large-distance behavior of correlation functions 
for all rj below the critical value 1/4, since large-distance 
fluctuations of the phase are still described by the Gaus- 
sian theory. Using that together with Eq. (|^5)) . we will 
now obtain the prediction for the spectrum of density 
ripples which is valid as long as only points with relative 
distances much larger than A2D contribute significantly 
to the integral in Eq. We will show below, that this 
regime is realized if 



> A 



2D- 



We introduce a dimensionless variable 



(49) 



(50) 



and use expression Eq. (j45p for all r. Using symmetries 
of the resulting integral, the expression for (|p(g)P) is 
written as 



2D 



ht 



l-r, 



F{v,y), (51) 



where i^(?7, y) is a dimensionless function defined by [12 

4 r°° r°° 

y Jo Jo 



{tx + yY + rlJ [tj. - vY + rl 



rl+rl 



(52) 



We find that the spectrum of density ripples remains 
self-similar in the course of expansion and the shape of 
the spectrum is a function of 77 only. Plots of Fijj^ y) 
for three different values of r] are shown in Fig. [51 and 
have a similar structure. Positions of maxima and min- 
ima are very well described by Eq. pg]) . where the up- 
per (lower) sign corresponds to the nth maximum (mini- 
mum). In Eq. (j23p typical distances which contribute to 
{\p{q)\^) near its maximum at y « tt can be estimated as 

Note however. 



^ y^ht/m, which leads to condition 
that self-similarity starts breaking down for sufficiently 
large y even when condition ([^5)1 is satisfied. 
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FIG. 6: Dependence of universal functions ^(77, y) on y — 
plotted for three different values of correlation exponents 
rj. Under condition (|49p functions -^(77, y) determine the self- 
similar shape of the spectrum of density ripples according to 
Eq. (I51|l . Curves from top to bottom correspond to ?7 = 0.25 
(solid, the Berezinskii-Kosterlitz-Thouless point), i] — 0.15 
(dashed), and 77 = 0.10 (dotted). 



Scaling of the magnitude of {\p{q)\'^) with time in the 
self-similar regime can be used to extract rj. For example, 
the integral of from zero to its first minimum 

scales with time as 



dq{\p{q)\')^t'/-'-\ 



(53) 



and the exponent changes considerably as 77 changes from 
to the critical value 1/4. 

For small 77, one can derive an expansion of ^(77, y) as 



where fi{y) can be evaluated analytically as 
/i(y) = 27rsin2|. 



(54) 



(55) 



The term f2{y) leads to a finite value of F{ri,y) at the 
first minimum. By including effects of 72(2/) and fz{y)i 
one can derive 

-J^^ ~ ^ (1-19^ + 0.38772) for 77 « 1, (56) 

which coincides with the direct numerical evaluation up 
to 2.5% for 77 = 0.25. 

For weakly interacting uniform 2D Bose gases at low 
temperatures, one can also obtain predictions which are 
not limited by Eq. (I49p . Under condition 



(57) 



an extension of Bogoliubov theory to 2D quasiconden- 
sates describes correlations at all distances [73|- Such a 
theory is valid up to temperatures of the order 



■log- 



n2D£,2D > 1; 



(58) 



and predicts the exponent (|47p . The correlation function 
is written as 



51 (r) 



7/;t(r, 0)7^(0,0) 



n2D 



exp 



27r7^ ^ / 



n2D£.2D > 1 



where the dimensionless function /2Z)(s) is defined by 

krlk 

Here, Jo{x) is the Besscl function, and Uk,Vk and Uk are 
defined by Eqs. (l28l) - ([30ll . 

We now consider a case of ®^Rb atoms with transverse 
confinement frequency uj± = 2tt x 3 kHz, and density 
77-2D ~ 84/zm~2. This yields the dimensionless interac- 
tion parameter g sa 0.13, healing length ~ 0.3 /im, 
and n2D^2D ~ perform a numerical integration 

of Eq. ((23)) for temperature T = 60 nK (which corre- 
sponds to correlation exponent 0.02) and various expan- 
sion times, and the results are shown in Fig. [T] 

Qualitatively, they look similar to the self-similar 
regime for all times, and one again obtains an oscillat- 
ing spectrum of density ripples with maxima and minima 
very well described by Eq. (|39p . In the weakly interacting 
regime the ratio of the first maximum to the first mini- 
mum for {\p{q)\^) is much larger than one, similar to the 
weakly interacting ID Bose gas. 



V. CONCLUSIONS 

To conclude, we calculated the evolution of the two- 
point density correlation function of an ultracold atomic 
Bose gas released from a tight transverse confinement. 

For ID gases in the weakly interacting regime, in a 
wide range of parameters given by Eq. (j33|) . we analyti- 
cally calculated the spectrum of density ripples (|p(g)P) . 
Our results are summarized in Eq. ([57)1 and Figs. [1] and 
[21 Our analytical theory is also applicable in the quasi- 
ID regime when ksT and /i are of the order of transverse 
confinement frequency ?iu;± . For expansion times smaller 
than 6.5mA|n/7i, we find that the spectrum of density 
ripples can have several maxima and minima, and their 
positions can be estimated using Eq. (|39p. While po- 
sitions of maxima and minima are essentially indepen- 
dent of the temperature, their amplitude exhibits strong 
temperature dependence. For ID quasicondensates, the 
density profile in external harmonic confinement depends 
weakly on the temperature, when the latter is of the or- 
der of the chemical potential [t^ • The density profile fol- 
lows the inverted parabola shape [zll, thus the bimodal 
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FIG. 7: (Color online) Normalized spectrum of density rip- 
ples (ip(fi')P) /{n2D^2D) for weakly interacting 2D quasicon- 
densate of *^Rb atoms with density n2D ~ 84/im~^, trans- 
verse confinement frequency lu± = 2n x 3 kHz, healing length 
~ 0.3 /im, and dimensionless interaction parameter g « 
0.13. Temperature is taken to be T = 60nK(A:Br ~ n), 
which corresponds to correlation exponent rj — 0.020. Vari- 
ous curves correspond to expansion times t — 12 ms (blue, 
solid), t = 8ms (green, dashed), and t — 4ms (red, dotted). 
The lines are guides to the eyes. 



density fitting cannot be used to measure temperatures 
reliably. We propose that our analytical result Eq. ([57|) 
can be used for thermometry of one-dimensional systems. 
Experimental investigation of this question is currently 
unde r wa y, and will be presented in a separate presenta- 
tion [l26| . 

For one-dimensional systems, we also discussed evo- 
lution of the density correlation function in real space, 
g2{x]t). For long expansion times we find that the cor- 
relation function g2{x;t) reaches the value of 2 at short 
distances and approaches the value 1 for distances larger 
than the correlation length, see, e.g.. Fig. [3] for Tonks- 
Girardeau regime. 

For 2D Bose gases with temperatures below the 
Bcrezinskii-Kostcrlitz-Thouless transition and suffi- 
ciently long expansion time, we showed that the spectrum 
of the density ripples evolves in a self-similar way. Our 
result for this case is given in Eq. (|?T|) and Fig. [HI with 
positions of maxima and minima determined by Eq. (j39p . 
The scaling of the overall magnitude can be used to ex- 
tract the correlation exponent 77, e.g., using Eq. ([55)) . 

For more complicated situations, e.g., multicomponent 
gases, relation (PT|) and its cross-correlation generaliza- 
tions can be used as a convenient experimental tool to 
characterize complex many-body states and their corre- 
lations. In addition, it can be used as an experimental 
tool to investigate non-equilibrium phenomena in low- 
dimensional gases. 
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APPENDIX A: TWO-PARTICLE DENSITY 
MATRIX OF A STRONGLY INTERACTING ID 
BOSE GAS 

In this appendix we will describe a Fredholm-type de- 
terminant representation for g{xi;x2',0) for < a:i < 
X2, which can be easily evaluated numerically. Due to 
Eqs. pS]) and this defines g{xi; X2; 0) for any values 
of xi and X2- Representations similar to the one devel- 
oped here can be obtained for any multi-point correlation 
function of bosonic fields in the strongly interacting limit. 

Mathematically, fermionization can be written as 



ip^{x) — exp 
ip{x) = exp 



in / dy%l)]{y)i}f{y) 



in / dyij]{y)ipf{y) 



^{x), (Al) 
V'/(x), (A2) 



where we introduced fermionic creation and annihilation 



operators '(/'^(a;) and ipf{x) 
commutation relations 



which have standard anti- 



[^{x)Jf{y)}^d{x~y), (A3) 

{^{x),^Pl{y)] - [i:fix)Jfiy)] = 0. (A4) 

For zero temperature, ground state for fermions corre- 
sponds to a filled Fermi sea, whereas at finite temper- 
ature one should use a thermal density matrix for non- 
interacting fermions. 

For convenience, we will introduce a fictitious under- 
lying lattice of spacing a ^ cci, 2:2, such that 



Xl 1 

— - = mi > 1, 
Z a 

X2 1 

— -- = TO2 > mi > 1, 
2 a 



(A5) 
(A6) 



where mi and m2 are large positive integer numbers. 
At the end of the calculation, we will take the limit 
a — > such that mio —f xi/2, m2a X2/2. On a lattice, 
fermionization rules [Eqs. (|Aip and (|A2p ] and commuta- 
tion relations (IA3I) and (IA4I1 are written as 



^^(*) - n [1 - 2^}(fc)V^/(fc)] (A7) 
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Vi(i) = n[l-2V'}(fc)Vi/(fc)] V-zW, (A8) 

k<i 

{^}W,V'/(fc)} =<5,fe, (A9) 
{V'}«,^}(fc)} - {V'/W,V'/(fc)} =0. (AlO) 
Using these relations, g{xi; X2', 0) can be written as 
£<(xi;x2;0) = ^?/'J-(TOi)?/'J-(-mi) 

X [] [l-2il,}{k)4'f{k)]iPf{m2)iPf{^m2)\, (All) 

where subset S equals 

S = [-7712 + 1, -mi - 1] U [mi + 1, ma - 1] . (A12) 
Expanding the parentheses, we obtain 



Since the structure of the matrix Aff""''^'' does not de- 
pend on n, summation over different n and sets ji, jn 
can be now r epre sented as a single Fredholm-type deter- 
minant Il29l . Il30l | 



Det [^y - 2aB,j] 



4a2 



(A18) 



where matrices Aij and -By of size 2(m2 — mi) x 2(m2 
mi) are defined by 



and 



A, = Diag{0,0,l,...,l}, 

Bij = G{si, tj), 



(A19) 
(A20) 



g{xi;x2; 0) = / ^(-2)"7/;t (m^V-Jl-mi) 

\ri=0 

X ^ijl)--1p}ijnHf{jn)--'tpf{jl) 

Jl<...<ir.Jfc6S 



X i^f{m2)i^f{-m2)) . (A13) 

For each n and set of ji, ■■■,jm expectation value of n-|-2 
creation and n + 2 annihilation operators can be written 
using Wick's th eore r n [l27l| as a determinant of (n + 2) x 
{n + 2) matrix [Hi, Il29l | 



(A14) 



where 



si = -mia, S2 = mifl, Si>2 = ji-2a, (A15) 
ti = m2a, t2 = -m2a, Uy2 = Ji-2a, (A16) 

and G{x, y) ~ G{x — ?/) is a Green's function of a free 
Fermi gas, which e.g. for zero temperature equals 



G{x) 



%xp[7M- = ^^^^^^- (A17) 

kf 27r TTX 



Si = —rriia, §2 — mia, (A21) 
Si — (— m2 + 7 — 2) a for 3 < 7 < 2 + m2 — mi, (A22) 
Si = (2mi — m2 + i — I) a for 
2 + m2 - m,i < i < 2(m2 - mi), (A23) 
ii = m2a, t2 = -m2a, ii>2 = s^. (A24) 



Expansion of the determinant of Aij — 2aBij using the 
rule for the determina nt o f the sum of two matrices (see, 
e.g., p. 221 of Ref. jl3l| ) generates the expansi on o f 
Eq. (jA13p , similar to a usual Fredholm determinant [l3Cl| . 
Indeed, only diagonal minors not including lines 1 and 2 
can be chosen from the matrix Aij . Complimentary minor 
of size (ri + 2) X (n + 2) from the matrix Bij is propor- 
tional to matrix M("+2) in Eq. (KT^i . and the summation 
over possible different sets of ji, is equivalent to a 
summation over different partitions of matrix Aij into 
diagonal minors. 

Since determinants are easy to evaluate numerically, 
one can now take the limit a numerically and eval- 
uate g{xi;x2',0) with any precision. 
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